REDUCTIONS AND DEVIATIONS FOR STOCHASTIC 
PARTIAL DIFFERENTIAL EQUATIONS UNDER FAST 
DYNAMICAL BOUNDARY CONDITIONS 

WEI WANG AND JINQIAO DUAN 

Abstract. In order to understand the impact of random influences at 
physical boundary on the evolution of multiscale systems, a stochastic par- 
tial differential equation model under a fast random dynamical boundary 
condition is investigated. The noises in the model and in the boundary 
condition are both additive. An effective equation is derived and justi- 
fied by reducing the random dynamical boundary condition to a simpler 
one. The effective system is still a stochastic partial differential equation. 
Furthermore, the quantitative comparison between the solution of the orig- 
inal stochastic system and the effective solution is provided by establishing 
normal deviations and large deviations principles. Namely, the normal de- 
viations are asymptotically characterized, while the rate and speed of the 
large deviations are estimated. 



1. Introduction 

The random fluctuations may have delicate impact in modeling, analyzing, 
simulating and predicting complex phenomena. The need to quantify un- 
certainties has been widely recognized in geophysical and climate dynamics, 
materials science, chemistry, biology and other areas [26j [HI [37]. Stochas- 
tic partial differential equations (SPDEs or stochastic PDEs) are appropriate 
mathematical models for various multiscale systems under random influences 



A stochastic partial differential equation usually contains noise in the equa- 
tion itself, i.e., the noise is acting on the system inside the physical medium 
[3H 121 El HO]- However, noise may affect a complex system not only inside 
the physical medium but also on the physical boundary. Randomness in such 
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boundary conditions are often due to various fast time scale environmental 
fluctuations. 

The usual boundary conditions, such as the Dirichlet or Neumann bound- 
ary conditions, do not contain time derivatives of the system state. On the 
contrary, dynamical boundary conditions contain time derivatives of the state. 

The boundary conditions may further contain random effects, as in some 
applications. For example, the environment surrounding a pipe fluid is usually 
subject to uncertain fluctuations, such as random vibration around a natural 
gas pipe or a waste water pipe. In a fluid laboratory, a wind tunnel or fluid pipe 
may be sitting on a flat foundation, which is also subject to random vibration. 
This noise affects the pipe fluid flow via boundary conditions such as a Dirichlet 
or Neumann boundary condition on a part of the pipe surface boundary, which 
is a static boundary condition perturbed by random fluctuations. The salinity 
flux on fluid inlet boundary of a gravity current (e.g., at the Strait of Gibraltar) 
has a fluctuating component and this leads to a random Neumann boundary 
condition OH]. Stochastic dynamical systems under such random boundary 
conditions have been studied recently in, for example, [TTj, [231 E21 ESI US, H]- 

In some other applications, the evolution of systems may also be subject 
to dynamical boundary conditions (containing time derivatives of the system 
state), under random perturbations. Such random dynamical boundary con- 
ditions arise in the modeling of, for example, the air-sea interactions on the 
ocean surface [33], heat transfer in a solid in contact with a fluid [27], chemical 
reactor theory [28] , as well as colloid and interface chemistry [39] . The random 
fluctuations on the boundary are usually much faster than the intrinsic time 
scale of these systems [15]. In these cases, the mathematical models are sto- 
chastic PDEs with fast-varying random dynamical boundary conditions. Such 
stochastic dynamical systems have been more recently investigated in, for ex- 
ample, pa U 033 M G2- 

Motivated by better analytical understanding of the above-mentioned mul- 
tiscale systems under fast scale random dynamical interactions on the physical 
boundary, as well as noisy forcing inside the physical medium, we consider 
a stochastic parabolic partial differential equation on a bounded domain D 
under fast varying random dynamical boundary condition on a part of the 
boundary dD. The fast time scale in the random dynamical boundary condi- 
tion is controlled by a small parameter e > and perturbed by a noise (white 
in time but correlated in space). Specifically, we study the following stochas- 
tic parabolic partial differential equation with a random dynamical boundary 
condition on a part of boundary I\ and a Dirichlet boundary condition on the 
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rest of boundary: 

du e = [Au e + f(u e ))dt + aidWi(t) , in D 
edu e = [—d u u e — u e ] dt + y/e &2 dW2(t) , on T\ 
u e = , on T 2 

where r\ and T 2 form the whole of the boundary dD of domain D, f(u) is 
some nonlinear term, and ax, 02 and e are constants. More details of this model 
will be presented in the next section. 

First, we derive an effective model for the above system as e — > . The 
effective model is still a stochastic partial differential equation, but with a 
simpler boundary condition (see Theorem 15.21) . Note that simple boundary 
conditions not only facilitate theoretical analysis but also are desirable for 
numerical simulations. To this end, we formulate the above stochastic sys- 
tem as an abstract stochastic evolution equation with non-Lipschitz nonlinear 
term having polynomial growth and fulfilling a suitable dissipativity condition. 
Since the nonlinear term is non-Lipschitz, we introduce a cut-off function and 
a stopping time to obtain a unique mild solution, which is also the unique weak 
solution for the system (see §3). Here we follow a semigroup approach which 
is also used in [B] to obtain a mild solution for a class of react ion- diffusion 
equations with multiplicative noise and non-Lipschitz reaction term but with 
deterministic static homogenous boundary condition. Next we present some 
useful a priori estimates for the weak solutions which yield the tightness of the 
distributions of the solutions (see §4). Then by a discussion in the variational 
form of the system we can pass the limit e — > to get the effective equation 
which is a stochastic partial differential equation but with a simpler boundary 
condition. And further we show that u e converges to the effective solution u 
of the limiting equation, in probability, in an appropriate function space. 

Then we determine the normalized deviation between the solution u e of the 
original system and the solution u of the effective system we obtained, it is 
proved that the normalized deviation 



converges, as e — > 0, to a process which solves a linear partial differential equa- 
tion with random coefficients, under a random boundary condition. Namely, 
this random boundary condition is a deterministic static boundary condition 
perturbed by a white noise; see Theorem 16.31 

Finally, we investigate the deviations u e — u of order e K , with < k < 1/2. 
In fact a large deviation result (see Theorem IT. 8j) is proved for 




[u e - u] 
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In this weak convergence approach we prove the Laplace principle which is 
equivalent to the large deviations principle in a complete separable metric 
space (i.e., Polish space). For background see [5] [T71 |2"T] . 

In [23] the authors have studied a system of reaction-diffusion equations 
with Lipschitz nonlinear term in open interval (—1,1) with a static boundary 
condition perturbed by a stationary random process which varies fast in time. 
The solution is represented through a Green function, then the limit e — > is 
passed in the space C([0, T] x (—1, 1)) and one gets a deterministic partial dif- 
ferential equation. The normal deviations of the solutions are then obtained in 
a weighted space with some assumptions on mixing properties of the random 
stationary process. 

This paper is organized as follows. The problem formulation and some 
preliminary results are presented in §2 . Section 3 is devoted to the derivation 
of some useful a priori estimates. The tightness of the laws or distributions 
for the solutions is proved in §4 and then the effective model is derived in §5 . 
The normalized deviations are studied in section §6 and the last section, §7, 
is devoted to a large deviations result. 



2. Problem formulation 

Let D be a bounded smooth domain in (1 < N < 3), with boundary 
To . Assume that r = Ti U r 2 , where T\ and r 2 are open subsets of To and 
Ii H r 2 = • We consider the following stochastic partial differential equation 
with a random dynamical boundary condition on Ti and a Dirichlet boundary 
condition on T 2 : 

du e = [Au e + f(u e )] dt + a 1 dW 1 (t) in D, (2.1) 

edu e = [ - d u u e - u e ] dt + ^a 2 dW 2 {t) on T 1 , (2.2) 

u e = on r 2 , (2.3) 

Ue (0) = u° in D , (2.4) 

7l u e (0) = 7iu° on ri (2.5) 

where e is real number with < e < 1 ; v = (n 1; n 2 , • • • , n N ) is the outer unit 
normal vector on T\ , and 

N 

d v = ^ nid Xi . 
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Wi(t,x) and W 2 (t,x) are mutually independent L 2 (D)-valued and L 2 (Ti)- 
valued Wiener processes, respectively, on a complete probability space (ft, J 7 , P) 
with a canonical filtration (J-"t)t>o • Moreover, 71 is a trace operator on Fx (see 
next paragraph), and o\ and a 2 are both constants. In fact, o\ := (J\Il 2 {D) 
and a 2 := cr 2 -Tk 2 (ri) , where Il 2 (d) and -7l 2 (Ti) are the identity operators on 
L 2 (D) and L 2 (Ti) respectively. Taking e = we formally have 

du = [Au + f(u)] dt + a 1 dW 1 (t) in D, (2.6) 

d u u + u = an Ti, (2.7) 

u = on T 2 , (2.8) 

M (0) = u° in D. (2.9) 

In the following we rewrite the equation in an abstract setting. Denote 
by H S (D) the Sobolev spaces W^D) , s > 0, with the usual norms; see for 
instance |30j for the definition. Note that H°(D) = L 2 (D) . We also define 
Hp 2 (/}) and Hq(L>) , s > 0, the spaces consisting of the function of H S (D) 
which vanish on the boundary T 2 and dD , respectively. In addition we denote 
by B s (ri) the Sobolev space H s -2(r!) , s > \ on the boundary T 1 . And let 
71 be the trace operator with respect to Ti which is continuous linear operator 
from H S (D) to B s (ri) for s > | . For more information on trace operators we 
refer to [38]. We also define spaces H~ S (D) , s > and B 1_s (ri), s > \ as the 
dual spaces of Hq(D) and B s (rx) respectively. 

We denote by (• , -)d , (• , -)ri the usual inner products in L 2 (D) and L 2 (Ti) , 
respectively. 

For our system we introduce the following functional spaces 

Xl = {(u,v) G Hj. a (D) x B 1 ^) : v = V~£Jiu} 

and X° = Lj, 2 (D) x L 2 (Fx) . Here L\-. 2 (D) is the space consisting of functions 
in L 2 (D) which vanish on T 2 . Define the norm and inner product on X° and 
X] respectively as 

\(u,v)\ 2 x0 = \u\ 2 L 2 (D) + \v\ 2 mri) , (u,v) e X°, 



and 



\l 2 (d) ' rlL 2 (ri) 
<(wi,t>i), (u 2 ,v 2 )) x0 = (ui,u 2 ) D + (vi,v 2 )r 1 

\(u,v)\ 2 x} = \u\ 2 Hl , D) + \v\l HTl) , (u,v) e Xl 



((tti,«i) , {u2,v 2 )) xi = (ui ,u 2 ) n i^ D) + (vi ,v 2 )bht 1 ) 
for [ui , Vi) G X° , i = 1,2. Here |m| h i ^ is taken as the equivalent norm as 
\u\nl (D) = |Vm| L 2 . 
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Remark 2.1. Notice that the space X\ depends on e but in the next section 
we just give some estimates for a fixed e > . And for passing the limit e — > 
we will consider in a fixed space instead of X\ . 

Now we define a boundary operator B on H 1 (D) as 

Bu = 7ictit + 7iu , ii G tffD) 

and a second order differential operator A = —A with homogenous Neumanna 
boundary condition. Then we introduce the operator A on D(A) = {(u,v) G 
X] : (-Am, Bu) G X } as 

Az = (An , -j=Bu) , z = (u ,v) G X\ . (2.10) 

Associate with the operator A we introduce the following bilinear form on X\ 

a(z, z) = (Az, z ) xQ = / VuVudx+ / ('jiu) ('jiu) dTi (2.11) 

J D JF! 

with z — (u , v) , z — (u, v) G X e x . Noticing that 

\lMWr)<C{T x )\u\^ D) , (2.12) 

there is some constants M > , a > and /3 G M such that 

a(z,z) < M\u\ n i {D) \u\ H i {D) (2.13) 

and 

a(z,z) > a\u\ 2 Hl(D) - p\u\ 2 L 2 [D) . 
Also by (12.121) then the following coercive property of a holds 

a(z,z) >a\z\ 2 xi - (3\z\ 2 x0 , z G X\ (2.14) 

for some constants a > and /3 G R. Then the linear operator —A generates 
a Co-semigroup, S(t), which is compact and analytic on X° ; see pp. 

Remark 2.2. Sometimes we also use the notation a(u,u) instead of a(z,z) 
for any z = (u,v) , z = (u, v) G X] . 



Now for the nonlinear term we make the following assumptions 

(F): / : R — > R is C 1 -continuous and there are positive constants a± , b\ 
such that 

f(u)u < —aiU A + bi , for u G R 
|/(u)| < aiH 3 + &i, forwGl 
f(u) < -a x u 2 + 6i , for ugR. 
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Define F(u) = f Q u f{y)dv. Then by the assumptions (F) there is some 
positive constant b such that 

F(u) < b and f(u) <b, VueR. (2.15) 

For the stochastic term we assume the following conditions. 

(£): Stochastic process W{t) = (W\(t) , W^if))* , is a Q- Wiener process 
on X° , defined on a filtered probability space (f2 , J 7 , T t , P) with co- 
variance operator Q = (Qi , Q 2 ) is trace class. Furthermore we assume 
Tr(AtQi) < oo. 

Remark 2.3. An example of such functions f is given by the following cubic 
polynomial 

f(u) = -au 3 + bu 2 + c 

with a > , b, cel. 

Remark 2.4. In (£), the technical condition on Qi is for the proof of regular 
properties of solution. As one example for such Qi , one can define W\ = 
y/ q(x)w(t) with positive function q E Hq(-D) and w a standard scalar Wiener 
process. Then, by the property of trace operator, the covariance operator of 
Wi , Qi = q satisfies Tr(A l l 2 q) = l^l^jip) < 00 ■ 

With the above notations system (I2.ip — (j2.5D can be written as the following 
abstract stochastic evolutionary equation 

dz e (t) = [ - Az e + H e (z e )] dt + Z dW(t) , z(0) = z° (2.16) 

where z e = (u e , \/e7i M e)*, H £ (z)(x) = (f(u(x)), 0))*, £ = (a u cr 2 )* and W(t) = 
(W\(t) , W 2 (t)y . The equation (12.161) can be further rewritten in the following 
mild sense 



7) 



z e (t) = S(t)z° + [ S(t-s)H e (z e (s))ds + [ S{t-s)ZdW(s). (2.1 
Jo Jo 

An adapted process z t is called a mild solution of (I2.16P if (12 . 1 7[) hold. For 

u e E L 2 (0,T;Hf 2 (D)) n C(0, T; L 2 (D)) , we call u e a weak solution of fl27lD - 

(1231) if for any t E [0,T) 

(u e {t),^{t)) D + e( llUe {t),^{t)) ri (2.18) 
= (u e (0),^(0)) 2J + e< 7l « e (0),V(0))r 1 + y (u e (s) , D ds 

+ e J ((7!tt e )(s), — ) Fi ds + j a(u e ,ip)ds + J (f(u e ),ip) D ds 

+ f ^,a 1 dW 1 (s)) D + v ^ [ (i>,o- 2 dW 2 (s)) ri , 
Jo Jo 
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for any if) G C 1 (0, T; C°°(D)). For more detail about solution of SPDEs we 
refer to [12]. 

We end this section by recalling the following two lemmas from |29j, which 
will be used in our later analysis. 

Lemma 2.5. Let Q be a bounded region iiiRxR". For any given functions 
h e and h in L P (Q) (1 < p < oo), if 

\h>e\LP(Q) < C , h e — > h in Q almost everywhere 

for some positive constant C , then h e — ^ h weakly in L P (Q) . 

Let X C 3^ C Z be three reflexive Banach spaces and X G y with compact 
and dense embedding. Define Banach space 



Lemma 2.6. If G is bounded in Q, then G is precompact in L 2 (0,T;y). 



In this section we prove the stochastic evolutionary equation (12. 16)) is well- 
posed and further derive a few useful a priori estimates on the solutions. Since 
the nonlinear term is non-Lipschitz, we apply the cut-off technique with a ran- 
dom stopping time. The same idea was used in [6j [7] for stochastic reaction- 
diffusion equations with local Lipschitz nonlinear terms and multiplicative 
noise. See also [3T] for systems on unbounded domain. 

Theorem 3.1. (Wellposedness) 

Assume that (F) and (S) hold. For any T > 0, let z° = (u°,v°) G X° be 
a \Tq, B(X )) -measurable random variable. Then system A2.16\) has a unique 
mild solution z e G L 2 (fl, C(0, T; X°) H L 2 (0,T; X})) , which is also a weak 
solution in the following sense 



G 



{h:hEL 2 (Q,T;X),— eL 2 (0,T; Z)} 



with norm 




3. Some a priori estimates 



(z e (t),<f>(t)) x o -(z e (0), <f>(0)) x0 



(3.1) 
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fort G [0,T) and G C 1 (0,T;X e 1 ). Moreover if z° is independent of W(t) 
with E|z°|^ < oo ; then there is positive constant Ct, which is independent of 
e, such that the following estimates hold: 

E\z e (t)\ 2 x0 + E f \z e (s)\ 2 xl ds <(l + E\z°\ 2 x0 )C T , fortG[0,T] (3.2) 
Jo 

and 

E{ sup \z e (t)\ x0 } < (1 + E|z°|^ + E / \z t (s)\ 2 xl ds)C T . (3.3) 
te[o,T] Jo 

Proof. For any integer n, we introduce the following cut-off function P n : IR + — > 
R + which is a smooth function satisfying P n (x) = 1 if x < n and P n (x) = if 
x > n + 1. Then nonlinear function P n (|z|)iJ e (z) is Lipschitz in both X° and 
X 1 where \z\ = \u\ + ^/e|7iit| for z = (u, y / e7i«). Now we have the following 
system with globally Lipschitz nonlinear term 

%(t) = S(t)z° + [ S(t-s)P n (\z?(s)\)H e (z?(s))ds+ [ S(t-s)ZdW(s). 
Jo Jo 

(3.4) 

Define a random stopping time r n (R) by 

T n {R) =inf{t>0: \ ~z n (t)\ > R] . 

Fix arbitrarily a positive number R < n and denote by xi the characteristic 
function of the set J. Consider the following integral equation for t < r n (R): 

rtAr n (R) 

z :(t) = s(t)z°+ s(t-s)P n (\z:(s)\)H e ( Z :(s))ds 

Jo 

rfAr"(fl) 

+ / S(t-s)ZdW(s). (3.5) 



Then by the Theorem 7.4 in [12], for any T > 0, the equation (13. 5p has a unique 
solution z n e G L 2 (Q,C(0,T;X°) n L^O^X 1 )). Moreover z™ is independent 
of n, foin>R and satisfies (I2.16P for t < r n (R). 

In the following we first derive some a priori estimates for z E in L 2 (Q, C(0, T; X°)n 
L 2 (0,T;X 1 )). Then we prove the wellposdness of the problem (12.161) . 

Applying the Ito formula to \z e \ x0 yields 

d\z e \ 2 x o + 2(Az e ,z e ) x °dt = 2{H e (z £ ),z £ ) x odt + 2(LdW(t),z e ) x o + \m\ 2 Q . 

Here II • II -q denotes the Hilbert-Schmidt norm of operator from Q l l 2 X Q to 
X° . Similar for || • || q 1 and || • |Lq 2 in the following. By the assumption (F) 
we have 

(H e (z e ),z e ) x o = (f(u e ),u e ) D < hMes(D) 
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where Mes(D) is the Lebesgue measure of the domain D. From the assump- 
tions (W) and (£), for any k G (1/2, 1), we deduce 

IIVl|2 _ 11^ 1 1 2 | li 112 

< a 2 trQ 1 + ajtrQ 2 . 
Now combining all the above analysis and ( 12.141) yields 

d\z € \ 2 x0 + a\z e \ xl dt < [b + a 2 trQ 1 + a 2 2 trQ 2 ]dt + 2{^dW{t),z,) xQ (3.6) 

where b = 2biMes(D). Integrating from to t and taking expectation on 
both sides of the above formula we have the estimate (I3.2p by the Gronwall 
inequality [12]. For a further estimate, we apply the Ito formula to \z, 
m > 1, and by calculation similar to get (13.21) we have 

j t v\ze\T + Efl^gr'kfti) < CiEi^ijr 2 + c 2 t 

for some positive constants C\ and C 2 - Then by the Gronwall inequality we 
have 

E|z e (*)|*S + E f \z t (s)\ 2 ™- 2 \z t (s)\ 2 xl ds < C T (1 + E|z°r*S). (3.7) 
Jo 

Integrating both sides of (I3.6p from to t, we have 

sup | z e (t)\ 2 x o < \z°\ 2 x0 + [aftrQ! + a 2 trQ 2 + b]T + 

0<t<T 



2m 
e\ X o, 



sup 

0<t<T 



(EW{s),z £ {s)) x ods +1. (3. 



By the Burkholder-Davis-Gundy inequality [13] and the assumption (S) we 
have 

-t 



E sup 

0<t<T 



(ZW(s),z e (s)) x ods 







< CE / \z e (s)\ 2 x0 \\E\\ 2 Q ds 

Jo 2 

< C'E [ \z e {s)\ 2 x0 ds 

Jo 

for some positive constants C and C. Then by (13.21) we obtain the estimate 
(CO) from (13~8D . 

Now we continue to prove the wellposedness of the system (12.161) . By the 
assumption of /, also satisfies the estimates (13. 2p and (13.31) which is inde- 
pendent of n and R. Then we have t(R) — > oo almost surely as R — > oo. For 
any T > define z e (t) = z™(t) for some n and R with n > R and r(i?) > T 
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almost surely. Thus P(r(i2) < T) = as R — > oo for any T > 0; for more 
details about proving global existence, see [3]. The uniqueness and continuity 
on initial value of z e follows from those of z™. Finally, by the stochastic Fubini 
theorem, i.e., Theorem 4.18 in [12], and the same discussion of |9J, we have 

The proof is complete. □ 

By Theorem 13.11 and the definition of z e we have the following corollary. 

Corollary 3.2. Under the same conditions as in Theorem \3.1l fort G [0, T], 
the following estimates hold: 

E(|w e (t)|i 2(D) + e| 7 i^(t)|i 2(ri) ) + 

j E(\u e (t)\n h{D) + ^iUe(t)g HVl) )ds < (1 + E|z°|^ )l7 t (3.9) 

and 

E{ sup [\uMh {D) + ^MhirJ} < (1 + E|^„)l7 t , (3.10) 
te[o,T] 

for some positive constant Ct , independent of e . 

Since the nonlinear term increases polynomially, in order to pass the limit 
e — > in system (I2.1l) - (l2.5p . we need a priori estimates for u t in the space 
Hp (£)) . In fact we have the following lemma. 



Lemma 3.3. Assume that (F) and (£) hold. Let z° = (u°,v°) e X° be a 
(JFq, B(X )) -measurable random variable with ~E\z \ xl < oo. Then for any T > 
0, the solution WIW z t = (u e , ^-nu € ) G L 2 (tt, L 2 (0, T; X 1 )nC(0, T; X )) and 
it satisfies the following estimates: 

EMt)|^ 2 ( D) +E^ \Au e (s)\l HD) ds < (l + E\z°\ 2 xl )C T (3.11) 

and 

E sup \u e (s)\ 2 nl (D) < (1 + E\z°\ 2 xl )C T (3.12) 
se[o,t] r2 

for any t G [0, T] and for some positive constant Ct, independent of e. 
Proof. Let V(z) = a(z, z). Applying the Ito formula to V(z e ), we have 

5S VW (3 - 13) 



-(A e z e ,A e z t ) Xo + (H(z e ),A e z e ) Xo + (TlW,A e z e ) 



X 



1 

+ 2 



l^lll 2 i + II^H^Qa 

/ -2 
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By definition of A e and (12.21) we have 
(H(z e ) , A e z e ) Xo = -(f(u e ),Au e ) D 

= (f'(u e )Vu e , Vu e ) D - (d u u e , f(u e )) Tl 
= (/'(« e )V« e , Vw e ) D + (eii e , /(u e ))n + (w e , /(w e ))ri 
-<Vi^ 2 ,/K))rx. (3.14) 
Also by the Ito formula 



(eit 



i e , /(u e )) ri = e J t ( F ( u t) > l)ri - ^-tr{D uu (F(u t ) , l) ri Qf (Qf )*) ■ 



Then we can have from (13.13j) - (13.14j) 
1 d 



2dt 



V(z e )-2e(F(u e ),l)r 1 



(3.15) 



|A^|xo - (ZW,AeZ e ) Xo + 



+ K 



-02 



+</ / (u 6 )Vu e , Vm £ ) d + (u e , /(u 6 )) ri - ^-tr{p uu {F{u e ) , I^QHqIT) ■ 

By assumptions (F) and (£), taking expectation on both sides of (13.151) . ap- 
plying the Cauchy inequality and noticing that 



u, 



|2 ^ (~i I |2 ill 2 

lL 2 (ri) S: ^e.\ u e\l?(p) + e Pe| H l (D) 



(3.16) 



for e > 0, we can have (13.111) by the Gronwall inequality and Corollary 13.21 by 
taking e > small enough. 

Integrating both sides of (13 . 1 5|) from to t, and using the Cauchy inequality, 
we conclude that 



1 



-E sup Mt)| H i 



te[o,T] 



< -E\u 

~ 2 



WAS 



\ L , ( D)ds+ I \u e {s)\ H x^ D) ds + 1 



for some positive constant C depending only on b, Qi, Q2, o\ and 02. Then 
we have (ETT^ by Corollary O- 

The proof is complete. □ 



4. Tightness of the distributions of solutions 

We intent to investigate the limit of the solution u e of stochastic system 
(I2.1l) - (l2.5p as e — > in the sense of distribution. For this purpose, in this 
section, we establish results on tightness of the distributions of the solutions. 
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Let /i e be the distribution of u e , which generates a Radon probability measure 
on the following metric space 

H := L 2 (0, T; L 2 (D)) n C(0, T\ H~ 1 {D)). 

Now we prove that the family of distributions {/i e } is tight in the space 7i. 

We apply the a priori estimates in the preceding section to obtain the tight- 
ness of {/i e }. 

First by the property of Wiener process, for some p G (0, |) 



E sup 

\t-r\<r 



\aiW(t) -aiW(r)| 



\t - r\P 



(4.1) 



Then, by estimate (13.21) and (I3.12p . for any given 5 > 0, there is a positive 

P{A 5 }> 1-5 



constant Cf. such that 



with 



\<j t W(t) - a t W{r)\ 



\t - r\P 



As = {ijj G Vt : sup Iwc^Ih 1 (d) < Cr> su P 

Q<t<T F ' 2 \t-r\<T 

and J \u € {t)\\ {D )dt<C 5 T }- 

For any ^ G Hq(D), by (J2HM23}, we have 

d{u e (s), ip) D = a(u e , <p)ds + (/(it e ), <p)r>ds + (aidW x {s), (p) D 
By the definition of a(z, z), we have 



L 2 (D) 



< C% 



(4.2) 



sup 

\t-r\<r 



a(u e , <p)ds 



< M 



\ U e(s)\ H i {D) ds 



1/2, 



vIh^d)^ (4-3) 



where M is defined in (12.1 3D . By the assumption (F) and the embedding 
HJ(D) into L 6 (L>) we have 



sup 

\t-r\<r 



{f(u £ ),ip) D ds 



< a 



■,1/2 



^Ifl-i^V^ ( 4 - 4 ) 



for some positive constant C" . Then if uj G As, by the definition of As and 

iD-flOD 

I M e I CP (0,T;H~ 1 (D))nL 2 (0,T;Hj, (D)) 

for some positive constant C(T, 5) . Define set Kg C 7i as 

K S = {u e H : |«e|c'>(0,T;ff- 1 p))nL 2 (0 ) T;Hi (£))) < C(T, 5)}. 
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Then by the compact embedding of C p (0, T; H~ 1 {D)) n L 2 (0, T; Hf 2 (L>)) into 
7i, Ks is compact in TC. And by the definition of As and above analysis we 
have 

P{u e G K s } >l-6. 
Then we proved the following result. 

Theorem 4.1. (Tightness of distributions of solutions) 

The family of distributions of the solutions, {fi e }, is tight in the space TC. 



5. Effective dynamics 

In this section we pass the limit of e — > in (12.1 8p and obtain the limit- 
ing system in a certain sense to be specified below. We always assume that 
E^ !^! < oo and E| 

^ 00 ^ n following. 
Since {fi e } is tight in the space TC (defined in the beginning of the last 
section), for any 5 > there is a compact set Ks C TC such that 

P{u e G K s } >l-6. 

Here Ks is chosen as a family of decreasing sets with respect to 6, i.e. Ks C Kg< 
for any 5 > 5' > 0. Moreover by the analysis of last section we can choose 
the set Ks with some positive constant Cy, depending on T and S, such that 
su Pte[o,T] I u t (t) | hi (D) < Ct for u e e K s . 

Then Prohorov theorem and Skorohod embedding theorem ([12]) assure that 
for any sequence {tj}j with tj — > as j — > oo, there exist subsequence { £?(£:)}, 
random elements {u* } cH,u* e TC, u°* G L 2 (D) and L 2 (D)-valued Wiener 

process Wf, L 2 (ri)-valued Wiener process W£ defined on a new probability 
space ¥*), such that 

C(u* , ) — C(u e .,.,) 

and 

u* , — > u* in TC as k — > oo, 

for almost all u; G fi*. Moreover u* solves system fl2.1|) - fl2.5p with W\ 
and W 2 are replaced by Wiener process Wf k and W^ k respectively defined on 
probability space (Q* , T* , P*) with same distribution as H^i and W2 for any /c. 
And u° is replaced by a random variable with C(u°) = C(u k *). And for 
P*- almost all u G Q*, \u° k * - u°*\ L 2 (D) and 

sup \W*(t) - W* k (t)\ L2(D) - 0, sup \W*(t) - W* 2k {t)\ LHTl) - 

0<t<T 0<t<T 
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for k — > oo. Now we will determine the limiting equation satisfied by u* and 
the limiting equation is independent of e. In fact we will prove that u* solves 
(12.6D — (12.9p with u° and W\ are replaced by u°* and W* respectively. 

We will pass the limit e — > in (12.181) for u* . For the nonlinear term f(u) 
increases polynomially, in order to pass the limit in f(u* ) we restrict u e in 

a bounded set in Hp AD). However this is impossible for stochastic process 
u* which con 

(fi*>|,P* 5 ) as 



u* which converges just in space 7i. For this define a new probability space 



and 



n; = W g n : «• w g 
.f; = {SnfiJ -.s e J 7 *} 



Denote by the expectation operator with respect to P|. It is clear that 
P(Q*\Q})<5. 

Since the distribution of u* is same as that of u e .,., , u* converges to u* 

in space L 2 (0, T; L 2 (D)) n C(0, T; H^p)) under the usual metric for lu G Cl* s . 
Here the usual metric in the intersection X f] Y, of two metric spaces (X, dx) 
and (Y, dy), is the metric d :— dx + <^y or equivalently d := max{(ix, ^y}- 

Now in the following we determine the limiting system satisfied by u* re- 
stricted on probability space (p,g,J-g,¥g). By the weak solution in the sense 
of (151) . for any fe^O.T; C°°p)) with ^(T) = we have 

" (u* e . w (t),iP) D dt-ej ( llU : m {t)^) Tl dt (5.1) 
a{u* e . {k) ^)dt+ J (f(u: m ),^) D dt + J (a 1 W* k (t)^) D dt + 

T 

(a 2 W; k (t)^) ri dt+(ul*,m)D + e( ll ul*, li m)r 1 , 

i 

on (ni,J7,PS) . 

We consider the terms in (15 .ip respectively. Since u*. converges weakly to 
u* in L 2 (0, T; Hf 2 p)) 

a(u* ,i/>)dt^> / a{u*,i))dt, e^O (5.2) 

for all a> G fix. By the definition of and the choice of \u* .(t)ln 1 m) is 
bounded uniformly in t G [0, T] and e G (0, 1]. By the embedding of Hp p) 

2iV 

into L W ^(D), \f(u* ) \l 2 (o, t-,l 2 (d)) is bounded. Then by Lemma [231 and 
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assumption (F), f(u* ) converges weakly to f(u*) in L 2 (0, T; L 2 (D)) which 

J (re) 

means for to £ tt* s 

T (fK j(k) ),4>) D dt^ j\f{u*)^) D dt, e^O. (5.3) 

By assumption (S) and the property of stochastic integral, see also Lemma 
3.1 of [25] 

Vi^r fc (s)^(*)Ms f \^Wl{s)^{s)) D ds (5.4) 



in PJ-probability and 



77 / (^m^Ms^O (5.5) 



in P^-probability. Moreover for u £ Q* s 



T 

(7i<, w (*)>^)ri<ft| 2 



v 

-> 0, e^O. (5.6) 

Then combining the above analysis in (I5.2l) - (l5.6p and by the density argument 
we could conclude 

fT 

(u*(t),ii>) D dt (5.7) 



a(u*,ip)dt+ / {f(u*),^) D dt+ / (criW*(t),^) r( dt + 
Jo Jo 

+ ( M °*,^(0)) D . 

Integrating by parts in (15 .7p we see that the limiting function w* satisfies 
the following system with deterministic boundary condition 

du* = [Ait* + f(u*)] dt + (7i dW*{t) in D, (5.8) 

d u u* + u* = on Ti, (5.9) 

u* = on T 2 , (5.10) 

u*(0) = u°* in D, (5.11) 

on the probability space (ft}, J 7 ^, P|), which has static boundary condition. For 
the above system fl5.8|) - fl5.11l) we can rewrite in the following abstract form 

du* = [-A u* + f(u*)]dt + o x dW?(t), u*(0) = u°* (5.12) 
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where the operator — A is the Laplace operator with the Robin boundary 
condition. The corresponding bilinear form is 

ao(vuv 2 ) = / VviVv 2 dx+ / (71^1) (71^2) dT t . 
Jd Jt! 



Then by Theorem 7.4 of |12] , and a similar analysis in the proof of The- 
orem 13.11 for any T > system fl5.14l) - fl5.17l) has a unique solution u* G 
L 2 (fl* s , L 2 (0, T; Hp 2 (D)) n C(0, T; L 2 (D))) in the sense of fETTl) . 
Then by the arbitrariness of the choice of 5, 

u* converges in P* — probability to u* (5.13) 

which solves 

du* = [Au* + f(u*)]dt + cti dW*(t) in D, (5.14) 

d u u* + u* = on Ti, (5.15) 

u* = on T 2 , (5.16) 

M *(0) = u°* in D, (5.17) 

on the probability space (fi*,jF*,P*). In order to obtain the convergence in 
probability of u e , we need the following lemma from [25J. 

Lemma 5.1. Let Z n be a sequence of random elements in a Polish space (X, d) 
equipped with Borel a-algebra. Then Z n converges in probability to an H-valued 
random element if and only if for every pair subsequences Z\ and Z m , there 
exists a subsequence Vk := (•Zj(fc), ^m(fe)) converging weakly to a random element 
v supported on the diagonal {(x, y) G X x X : x = y}. 

Here we call an X-valued random variable X n converges weakly to X if 

Ef(X n ) = [ f(x)F n (dx) - / f(x)P(dx) = Ef(X) 

with P n = C(X n ) and P = C(X). Notice that convergence in probability 
implies weakly convergence, see [2J IT5] . Then by the uniqueness property of 
solution for equations fl2.6l) - fl2.9|) which can be proved similarly by that for 
equations fl5.8p - ()5.1ip . we can formulate the main result in this section by 
Lemma 15.11 

Theorem 5.2. (Effective system) 

Assume that conditions (F) and (S) are satisfied. Let u° G L 2 (D) be a 
(JF , B(L 2 (D)))- measurable random variable, which is independent ofW{t), 
with Ej-u !^! < 00 and E[|m |^ 2 ( D ) + l7i M °ll2( ri )] < 00. Then for any 

T > 0, the solution u t of the stochastic system A2.1\) - I[27d)) converges to u, 
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which is the solution of the effective limiting system ^Uj\) - $2J$) , in probability 
in space H: 

U = L 2 (0, T; L 2 (D)) n C(0, T; H'^D)). 

6. Normal deviations 

We have proved u e approaches u in probability in space H, namely, the 
difference u e — u tends to in probability in space Ti as e — > 0. In this section 
we consider the order of u e — u in e as e — > 0, that is, the normal deviations 
of w e away from the effective solution u. We prove that the order is \ and the 
normalized difference ^(u e — u) converges in an appropriate function space. 
In this section we assume the assumption in Theorem 15.21 

Denote by v e = -^(u t — u) . Then we have the following initial boundary 
value problem for v e 

v e = Au e + i=(/(tt 6 ) -/(«)), in D (6.1) 

d u v e + v e = -\feu e + a 2 W 2 , on Ti (6.2) 
v e = 0, on T 2 (6.3) 
v e (0) = 0. (6.4) 

As e — * we expect v e converges in some sense to the solution v of the 
following linear system 

dv = [Av + f'{u)v]dt, in D (6.5) 

d v v + v = 0-2W2, on Ti (6.6) 

v = 0, on T 2 (6.7) 

v(0) = 0. (6.8) 

Note that the limiting system f l6.5p -( |6T8j) contains a static boundary with ran- 
dom force (but not dynamical). For the wellposedness of the above two systems 
we follow the approach of [TTJ [13]; see also [32]. Let AT be a linear bounded 
operator from L 2 (I\) into L 2 (D) defined as the solution of following problem 

ry — Ay = in D, d u y + y = g on Ti (6.9) 

with rei such that ( 16. 9j) has a unique solution y = Afg for any g E L 2 (Ti). 
Here Af is called the Nenumann mapping. For our problem let g = g e = 
— y / e , u e + 01W1 and go = o%Wi. However it is easy to see that g e and go are 
not in the space L 2 (Ti). Fortunately, we can extend Af to a bounded linear 
operator from i^Ti) to D £ A which is the domain of the operator (rl — A) 6 with 
< e < f + |, Q > — §, see [30] or example 3.1 in [32]. Here A is a second order 
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differential operator defined on {u G H 2 (D), d v u = 0} with Au = Au. Denote 
by S(t) the Cq semigroup generated by the linear operator A. Then we can 
write the solution of (I6.ip - fl6.4l) and (I6.5p - (I6.8I) respectively in the following 
mild sense 



v e (t) = - [ AS(t-r)j\f(^iV e )(r)dr + a 2 f AS(t - r)j\fdW 2 (r) + 
Jo Jo 

S(t-r)[f(u £ (r))-f(u(r))]dr + 



AS(t-r)j\f(-yiu e )(r)dr (6.10) 



o 



and 



v(t) = -[ AS(t-r)j\f{^ 1 v)(r)dr + a 2 [ AS(t - r)NdW 2 (r) + 
Jo Jo 



S{t — r) f (u(r))v(r)dr. 



o 



By the Example 3.1 of [32] for any T > 0, there exist functions v e and v , 
both in C(0,T; L 2 (n, L 2 (D))), which are unique mild solutions of fl6TT]) - fl6T4l) 
and f l6.5l) -( l6T8l) . respectively. 

Remark 6.1. For a special one-dimensional case on the domain D = (—1, 1), 
the solution is proved earlier [23] in a special weighted space C a C C a ([0, T] x 
(—1,1)) . Here C a ([0, T] x (— 1, 1)) consisting of all continuous functions u(t,x) , 
t G [0, T] , \x\ < 1, such that \iYn x _ r ± 1 a(x)u{t ) x) = uniformly in t G [0, T] . 
Here weighted function a(x) can be chosen as (1 — x 2 ) a , < a < 1 . Then 
C a C C a ([0,T] x (—1,1)) consisting ofu(t,x) such that h(t,x) = f Q u(s, x) ds 
has uniform limit in t G [0, T] when x — > 1 and when x — > — 1 . 



Let v e be the distributions of v t in the space L 2 (0, T; L 2 (D)). For our purpose 
in the following, we prove the tightness of z/ e . First we should derive a further 
a priori estimate for v e . As pointed out in [32J the Ito formula cannot be used 
for the Lyapunov function V(x) = \x\ p , p > 0. We treat v e in the mild sense 
(I6.10p . For any 5 > 0, we still consider u> G As which defined by (14. 2p . 

Lemma 6.2. Let u° G L 2 (D) be a (To,B(L 2 (D)))- measurable random vari- 
able, which is independent ofW{t), with EIm ^ ^ < oo and E[|m°|^ 2( . d ^ + 

l7i M °ll 2 (ri)] < 00 ■ Then for any T > 0, there exist a X > and a positive 
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constant Ct such that 

E [XA>e|£2 (0iT . H A 2 (D))] < c t- 

Here xa s {^) = 1 f or w G and xa s {^) — otherwise. 

Proof. By the similar estimates in the proof for Proposition 2.2 of |32j and 
estimates (13. 7p with m = p , noticing that the initial value is zero, we have 

E \.XA s \v £ \ q LHOtT . u ^ m ] <C T 

for some q > 2. Then by the Holder inequality we have the result. □ 
Furthermore multiplying if G Cq°(D) to both sides of (16. ip yields 

/ (Ve(t),<p) D dt 

Jo 

= - [ (Vv t (t),Vtp) D dt+±: f ((f(u e (t))-f«t))),Cp) D dt. 

Jo ' V e Jo 

Then by Lemma 16.21 we deduce 

nXA 6 \vMl H o,T^- H D))} < Ct (6.11) 

for some positive constant Ct- Then by the Chebyshev inequality and Lemma 
12.61 for any 8 > there is a compact subset Ng C L 2 (0, T; L 2 (D)) such that 

P{v e G N 5 } > 1 - 5. 

That is, the probability measure sequence {u e } is tight in space L 2 (0, T; L 2 (D)). 
Then Prohorov theorem and Skorohod embedding theorem Q12J) assure that 
for any sequence with 6j — > as j — > oo, there exist subsequence {ej(fe)}, 
random elements {v tj(k) } C L 2 (0, T; L 2 (D)), v G L 2 (0, T; L 2 (D)), L 2 (Ti)- 

valued Wiener process W 2 k defined on a new probability space (fl, JF, P), such 
that 

and 

v ej r k) — > v in L 2 (0, T; L 2 (D)) as k — > oo, 

for almost all G fi. u e , fe . solves (I6.ip - (l6.4p with W2 replaced by W^fc- And 
for almost all u> G f2 

sup |W 2 jfc - W 2 |l2(£») 0, k ^ oo. 

0<t<T 

Moreover tL.,„ = , 1 (■»,, — u) for some random elements u e . rus , u & H with 
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In order to pass limit e — > in /', by the same approach of §3 we define the 
following new probability space (^,5, J~5, IP5) as 

U 5 = {u' eU: v e (cu') G iV" 5 }, 

and 

= Pi(gn^) 

Now we restrict the system ( 16. ID — (16.41) on the probability space (0,5, T5, Pj). 
By the definition of fia and the discussion in $5] for almost all u G Qs 

u ej(h) (t, x) — > u(t, x), almost everywhere in [0,T]xD. (6-12) 
And for u G fig, ^e j(A .) converges to v almost surely on [0,T] x D, 
v ej , k) —> v weakly in L 2 (0, T; Hp 2 (D)), as — > 00 

and 

u e ^„. — > v strongly in L 2 (0, T; L 2 (D)), as k — > 00. 
Taking ^ G C7 CO (0,T; l7°°(D)) as the testing function for fl£H) yields 

T 

<w ej(fe) (s),^(s)) D rfs (6.13) 

T /-T 

(s), A^j(s)) D ds- / (W(s),n ej(i) (s)) r ds + 

1 

7^ 



</ 



[f(ue m (s) - f(u(s)),ip) D ds- / (v^ (fc) (s),^(s)) ri ds + 



./O 

<U fi . w (s),'0(s)) ri ds+ / <CT 2 H / 2fc(s),^(s)) ri C/s. 



'0 </0 

We pass the limit e — > in (16.131) . Notice that 



V^m / (^ (fc) (s),V(s)> ri ^ 

Jo 

< v/^iW / l^ j(fc) (s)|L2(r 1 )|V'(s)|L2(r 1 )^ 
Jo 

7o 2 ./o 2 



0, e — > 0, for c<j G fL>. 
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By the assumption (F) and (16.121) 
1 



U'( u ^ ik) (s)) - f{u(s))) = f(u em (s))v, 



with u 6j(k) — > u almost surely on [0, T] x D. Then by Lemma 12.51 f'(u €j(k) ) 
converges weakly to f'(u), for uj G And by the choice of v e ., k) , which 
converges strongly to v in L 2 (0,T; L 2 (D)), we have f'(u e . (k) (sj)v e . {k) converges 
weakly to f'{u)v in L 2 (0,T; L 2 (D)) which yields 

1 /" T ' T 
V^iW ./o 
Also by Lemma 3.1 of [25J 



oo 



(f(u e . (k) (s)) - f(u{s)),tp) D ds -> / (f'(u)v,ip) D ds for G ^ 



((J2W 2 k(s),ip(s)) ri ds (a 2 W 2 (s),ip(s)) ri ds, k 



in P^-probability. 

Then combining all the above analysis for the terms in (I6.13p . we can pass 
the limit e — > in (16. 13[) and conclude that 

T 

(6.14) 

Af(s),^(s)) ds+ / (f'(u)v,ip(s)) D ds 
Jo 

- I (v(s),if>(s)) ds+ I (<r 2 W 2 (s),i>(s)) ds, 
Jo Jo 

which is the variational form of (16.5p — (16.8p . Notice that we have proved the 
wellposedness of (I6.5p - (I6.8I) . Then by the arbitrariness of 5 and the same 
discussion in the proof of Theorem 15. 2} we have the following result on normal 
deviations. 

Theorem 6.3. (Normal deviations principle) 

Assume that the conditions (F) and (S ; ) are satisfied. Let u° G L 2 (D) be a 
(J-'o, B(L 2 (D)))- measurable random variable, which is independent ofW(t), 
with Elw !^! ^ < oo and E[\u°\^ 2 ^ + I7i w °ll2( ri )] < oo. Let u € and u be 

the unique weak solutions of $2. fixity and Il2.6]) - I[27ty) . respectively. Then 
-^(u e —u) converges in probability to a stochastic process v, which is the solution 

of the linear random system Ii6.5\) - I[67^) . in the space L 2 (0,T; L 2 (D)). 



SPDES UNDER FAST DYNAMICAL BOUNDARY CONDITIONS 



23 



7. Large deviations 

In §0 Theorem 15.21 we have proved that u t —>■ u as e — ► 0. We have also 
obtained convergence result of the normal deviations of order in £j6j Theorem 
16.31 which implies the normal deviations of order e K tend to for < k < |. 
In this section we consider the logarithmic asymptotics of the deviations of 
order e K , < k < |, in probability. That is, we consider the deviations of 
v* = e~ K (u e — u) which satisfies 

v: = Av: + e- K (f(u e )-f(u)), in D (7.1) 

d u v« + v« = -e 1 - K u e + e^~ K a 2 W 2 , on T 1 (7.2) 

v K e = 0, on T 2 (7.3) 

<(0) = 0. (7.4) 

We intend to prove that the family {v* : e > 0} satisfies the large deviations 
principle in L 2 (0, T; L 2 (D)). We follow the results on large deviations in [5] 
for Polish space valued random elements; see also [36] for large deviations of 
two-dimensional stochastic Navier-Stokes equations. 

Let Ho C H be Hilbert spaces with norm \-\h and | • \h respectively. Assume 
that the embedding of H in H is Hilbert-Schmidt. Define A the class of H - 
valued .^-predictable process w satisfying J Q T \w(s)\ 2 Ho < oo a.s. For M > 
let 

S M = {weL 2 (0,T;H ) : / \w(s)\ 2 Ho ds < M} 

Jo 

which is a Polish space (i.e., complete separable metric space) endowed with 
the weak topology. Define Am = { w £ A : w G 5m, a.s.}. 

Let .E be a Polish space and g e : C(0, T;H) —> E be a measurable map. 
Let 1/ be an iJ-valued Wiener process. Define Y e = g e (V(-)). We consider 
the large deviation principle for Y e as e — > 0. Since E is a Polish space, the 
Laplace principle and the large deviation principle are equivalent [36] . 

Definition 7.1. A function I mapping E to [0, oo] is called a rate function if 
it is lower semicontinuous . A rate function I is called a good rate function if 
for each M < oo, the level set {y e E : I(y) < M} is compact in E. 

Recall that a family {Y e : e > 0} of i?-valued random elements is said to 
satisfy the large deviations principle (LDP) with speed a(e) — » oo, as e — > 
and rate function J if (see [22| [12] ) 

(1) For any S, 7 > and y £ E, there exists eo > such that for any 
e€(0,e ) 

P - < 5} > exp{-a(e)(/(y) + 7)}. 
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(2) For any r, 5, 7 > 0, there exists eo > such that for any e G (0, eo) 
P - $(r)| B > <5} < exp{-a(e)(r - 7)} 

where $(r) = {y E E : I(y) < r}. 

It is well known that the large deviations principle and the following Laplace 
principle are equivalent in Polish space. 

Definition 7.2. Let I be a rate function on metric space E. A family {Y £ : 
e > 0} of E-valued random elements is said to satisfy the Laplace principle on 
E with rate function I and speed a(e) — > 00, e — > 0, if for each real valued, 
bounded and continuous function h defined on E, 

Km logE {exp [-a(e)h(Y*))} = - inf {%) + I(y)} . 

For our purpose we introduce the assumption 

(H): There exists a measurable map g° : C(0,T; H) — > E such that 

(1) Let {w e : e > 0} C Am for some M > 0. Let w e converges in 

distribution to w. Then g € I V(-) H — -)= f n w e (s)ds I converges 

\ y a ( e ) u / 

in distribution to g°(J Q w(s)ds). 

(2) For every M < 00, the set Km = {9°(J w ( s )ds) : w G <Sj\/} is a 
compact subset of E. 

For each g & E, define 



1 



T 



I(y) = inf 1 - / ds ^ . (7.5) 

{ W ^{Q,T;H ): y=g»{J oW (s)ds)} [2 J J 

Then we have the following theorem 

Theorem 7.3. Let Y e = g e (V(-)). If g e satisfies the assumption (H) ; then the 
family {Y e : e > 0} satisfies the Laplace principle in E with rate function I 
given by ^7.5\ ) and speed a(e). 

The proof of the above theorem is similar to that of the proof of Theorem 
4.4 in [5] which is for the speed a(e) = e _1 . We omit it here. 

In the following we apply the above result to the system (I7.1l) -( l7~4l) . In 

this case H = L 2 (r0, H = Q^H, E = L 2 (0,T; L 2 (D)), V(-) = W 2 (-) and 
a(e) = e 2 * 1 " 1 , < k < |. Since Q2 is a trace class operator, the embedding 
of H in H is Hilbert-Schmidt. By the analysis of Section [6] there exists a 
Borel measurable function g 6 : C(0,T;H) — > E, such that v* = g t (W 2 ). We 
intend to verify the assumption (H) for g e . In fact four lemmas are proved 
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to complete the verification. Let g e be defined as above. For any w G Am , 
< M < oo , denote g e (w 2 (-) + e~^ +K J Q w(s)ds) by . 

Lemma 7.4. is the unique weak solution of the following stochastic system: 

b K t = Av: + e-*(f(u e )-f(u)), in D (7.6) 

d v v K e +v K e = -e 1 - K u e + a 2 w + e^- K a 2 W 2 , on T 1 (7.7) 
v K t = 0, on T 2 (7.8) 
v*(0) = (7.9) 

in L 2 (n,C(0,T;L 2 (D))nL 2 (0,T;H£ 2 (D))). Here A is chosen in Lemma\EM 

Proof. This result follows from a Girsanov argument. In fact let W 2 (-) = 
W 2 {-) + e~5 +K J o w{s)ds. Then W 2 is Wiener process with covariation Q under 
the probability P w which satisfies 

dP w = exp je^ +K jT w(s)dW 2 (s) -^e- 1+2K \w(s)\ 2 H ds^ dP. 

Then a similar analysis in Section [6] yields the result. □ 

Lemma 7.5. Lei w G L 2 (0,T; H) . Then the following stochastic system 

p w = Ap w + f'(u)p w , in D (7.10) 

duPw + Pw = cr 2 w, on r x (7.11) 

p w = 0, on T 2 (7.12) 

p w (0) = 0. (7.13) 

has a unique weak solution p w G C(0,T; L 2 (D)) flL 2 (0,T; Hp 2 (£>)). 

Proof. This is a classical result of nonhomogeneous boundary problem [30J . □ 



We now define the function g° as follows: g°(h) := p w if h — j Q w(s)ds for 
some w G £ 2 (0, T; if), otherwise g°(/i) = 0. By the same discussion in Section 
[6] for the normal deviations we conclude that 

Lemma 7.6. Let {w e } C Sm converge in distribution to w, as a Su-valued 
random variable. Then g e (w 2 {-) + e K ~5 J w t (s)ds^ converges in distribution 
to g° (J Q w(s)ds) in E. 

Lemma 7.7. Let < M < oo be fixed. Then the set K M = {g°(J w(s)ds) : 
w G Sm} is compact in E. 
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Proof. By the definition of Sm, for any sequence {w n } C Sm, there is a sub- 
sequence w n (relabelled by n ) and w G Sm such that w n weakly converges to 
w as n — > oo. Then it is enough to prove that p w n converges to p w in E. Let 
@n = Pw" - Pw, then 

AQ n + f'(u)Q n , m D 
a 2 (w n — w) , on Ti 
, on r 2 
0. 

A simple energy estimate and the fact that embedding of Hq in H is Hilbert- 
Schmidt yield 

SUp \Qn(t)\ 2 L 2 {D) + [ |e n (s)|^i , D) ds 
0<t<T Jo 2 

< C(T,b,a 2 ) [ T \w n (s)-w(s)\ 2 L2(ri) ds 
Jo 

— > 0, as n — > oo. 

This completes the proof. □ 
By the Lemma [7.4H7.7I and Theorem 17.31 we can draw the following result. 
Theorem 7.8. (Large deviations principle) 

Assume that the conditions (F) and (S') are satisfied. Let u° G L 2 (D) be a 
(JFq, B(L 2 (D))) -measurable random variable, which is independent of(Wi(t), 
with Ej-u !^! < oo and E [Im !^^) + |7iM°|^2( ri )] < oo. Let u e and u be the 

unique weak solutions of h2.1\) - UZlfy and li5.8\) - l[5Jl\) . respectively. Then for 
any < k < | ; e~ K {u e — u) satisfies large deviations principle with good rate 
function /(•) given by \7.5\) and speed e 2K ~ l in L 2 (0,T; L 2 (D)). 



e n = 

d u e n + e n = 

e n = 

e„(o) = 
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